Dynamical properties of model communication networks 
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We study the dynamical properties of a collection of models for communication processes, char- 
acterized by a single parameter £ representing the relation between information load of the nodes 
and its ability to deliver this information. The critical transition to congestion reported so far oc- 
curs only for f = 1. This case is well analyzed for different network topologies. We focus of the 
properties of the order parameter, the susceptibility and the time correlations when approaching 
the critical point. For £ < 1 no transition to congestion is observed but it remains a cross-over from 
a low-density to a high-density state. For £ > 1 the transition to congestion is discontinuous and 
congestion nuclei arise. 

PACS numbers: 89.20.Ff,05.70. Jk,64.60.-i 



I. INTRODUCTION 

The interaction between the elements of social, tech- 
nological, biological, chemical and physical systems usu- 
ally defines complex networks. The study of topological 
properties in such networks has recently generated a lot 
of interest among the scientific community jjj [| Q, || . 
Part of this interest comes from the attempt to under- 
stand the behavior of technology based communication 
networks such as the Internet ||, the World Wide Web 
M, e-mail networks Q or phone call networks JlO| . 
The study of communication processes is also of inter- 
est in other fields, notably in the design of organizations 
JTT| , [l2[ . It is estimated that more than one-half of the 
U.S. work force is dedicated to information processing, 
rather than to make or sell things in the narrow sense 

Tools taken from statistical mechanics are used to un- 
derstand not only the topological properties of these com- 
munication networks, but also their dynamical proper- 
ties. Particularly interesting is the phenomenon of con- 
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gestion. In has been observed, both in real networks 
and in model communication networks, |14 
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that the system only behaves efficiently when the amount 
of information handled is small enough. The network 
collapses above a certain threshold and some informa- 
tion is accumulated and remains undelivered over large 
time periods — or it is simply lost. This transition from 
a free to a congested regime is indeed a phase transition 
and could be related to the appearance of the 1// noise 
observed in Internet flow data |l^] . 

Understanding congestion is also important from an 
engineering point of view. For instance, in October 
of 1986, during the first Internet collapse reported in 
the literature, the speed of the connection between the 
Lawrence Berkeley Laboratory and the UC at Berkeley — 
separated by 200 meters — dropped by a factor of 100 
p3| . Although the problem of congestion, and in par- 
ticular its prevention and control Ijjfl, has been studied 



because of its implications in digital communications, a 
deep understanding of the physics of congestion for gen- 
eral communication processes and beyond particular pro- 
tocols is still lacking. 

A general collection of models that captures the es- 
sential features of communication processes has been re- 
cently proposed |L6]]. In these models, agents — nodes — 
are organized in a hierarchical network — a Cayley tree — 
and interchange information packets. Each agent has a 
certain capability that decreases as the number of pack- 
ets to deliver increases. When the capability is inversely 
proportional to the number of accumulated packets, a 
continuous phase transition was found between the free 
and the congested phases. This transition was character- 
ized by means of an order parameter. 

The aim of this paper is to study the collection of mod- 
els mentioned above, since only hierarchical lattices and a 
very particular congestion behavior — capability inversely 
proportional to the number of accumulated packets — 
were considered. An extension to consider the existence 
of cost to maintain communication channels has already 
been done | pf]| . First, the network model is extended 
beyond purely hierarchical Cayley trees. In particular, 
1-dimensional and 2-dimensional lattices are considered. 
It should be noted that simple models such as Cayley 
trees and regular lattices can capture the main char- 
acteristics of the dynamics of information flow in com- 
plex networks. Hierarchical trees have been considered 
to model the TCP/IP protocol in the Internet @ 
On the other hand, computer based communication net- 
works have been described by placing routers and servers 
in square lattices JlJ, |i~7|| . 

Moreover it is shown that, independently of the topol- 
ogy, the collection of models can be split into three 
groups according to how the network collapses. In the 
first group, agents deliver more packets as they are more 
congested — although their capability, as defined in |L6[ , 
decreases — , and the network never collapses. In the sec- 
ond group, agents deliver always the same number of 
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(a) (b) (c) 

FIG. 1: Network topologies studied, (a) ID lattice with L — 
5; (b) 2D lattice with L = 5 and S = 25; and (c) (2, 3) Cayley 
tree with branching z = 2 and m = 3 levels or generations. 

packets independently of their load — number of packets 
to deliver — ; this behavior leads to a continuous phase 
transition as reported for hierarchical lattices in |L6| . Fi- 
nally, when agents deliver fewer packets as their loads 
increase, the transition to the congested phase is dis- 
continuous and the network collapses in an inhomoge- 
neous way giving rise to congestion nuclei. To character- 
ize these behaviors, the order parameter defined in ]l6| is 
used, but also the power spectrum of the fluctuations and 
a susceptibility-like function. Most of the current effort 
is focused on the quantitative study of the continuous 
transition and the critical behavior associated to such a 
transition. 

The paper is organized as follows. Section 2 presents 
the collection of models and describes some details con- 
cerning the extension beyond hierarchical lattices, includ- 
ing ID and 2D lattices. In section 3, these models are 
studied numerically and analytically. The scaling of the 
critical point with the size of the system, the behavior 
of the order parameter, and the divergence of the char- 
acteristic time at the critical point are studied in detail. 
Finally, section 4 includes the discussion of the results 
and the conclusions. 



II. THE MODEL 

The model considers three basic components: (i) the 
physical support for the communication process — agents 
and communication channels — , (ii) the discrete informa- 
tion packets that are interchanged, and (iii) the limited 
capability of the agents to handle such packets. 

The communication network is mapped onto a graph 
where nodes mimic the communicating agents (for in- 
stance, employees in a company, routers and servers in 
a computer network, etc.) and the links between them 
represent communication lines. In particular, the three 
different topologies depicted in Fig. |l| are considered: ID 
lattices of length L, 2D square lattices of side length 
L and size S = L 2 , and hierarchical Cayley trees with 
branching z and a number of generations or levels m, 
hereafter denoted (z,m). 

The dynamics of the model is as follows. At each time 
step t, an information packet is created at every node 
with probability p. Therefore p is the control parame- 
ter: small values of p correspond to low density of pack- 



ets and high values of p correspond to high density of 
packets. When a new packet is created, a destination 
node, different from the origin node, is chosen randomly 
in the network. Thus, during the following time steps 
t + 1, t + 2, ... , t + T, the packet travels toward its des- 
tination. Once the packet reaches the destination node, 
it is delivered and disappears from the network. 

The time that a packet remains in the network is re- 
lated not only to the distance between the source and 
the target nodes, but also to the amount of packets in 
its path. In particular, at each time step, all the packets 
move from their current position, i, to the next node in 
their path, j, with a probability q^. This probability q^ 
is called the quality of the channel between i and j, and 
was defined in |[16| as 



q%j = y/kikj , (1) 

where ki represents the capability of agent i at each time 
step. The quality of a channel is, thus, the geometric 
average of the capabilities of the two nodes involved, so 
that when one of the agents has capability 0, the channel 
is disabled. High qualities (qij sa 1) imply that packets 
move easily while low qualities (qij w 0) imply that it 
takes a long time for a packet to jump from one node to 
the next. The general equation proposed for ki is: 

h = f(rn) (2) 

where n, is the total number of packets currently at 
node i. The function f(n) determines how the capabil- 
ity evolves when the number of packets at a given node 
changes, 

, , , _ / 1 for n = , , 

fW ~ \ n"« for n= 1,2,3,... (6) 

with £ > 0. Equation (|^) defines a complete collection of 
models with agents that behave differently depending on 
the exponent £. 

The average number of packets delivered during one 
time step by a node i to another node j is proportional 
to ni/(n\ Uj ). Assuming that n,; oc nj the former 

expression reads . The proportionality is exact both 
in the hierarchical lattice and in large enough ID and 2D 
lattices, where adjacent nodes are statistically equivalent. 

Therefore, for £ < 1 the number of delivered packets 
increases with the number of accumulated packets. For 
£ > 1 the number of delivered packets decreases as the 
number of accumulated packets increases. Finally, for the 
particular case £ = 1, the number of delivered packets is 
independent of the number of accumulated packets. Note 
that this particular case is consistent with simple models 
of queues Q . 

The last point that needs to be explained to com- 
pletely define the current model is the routing algorithm 
or, in other words, the set of rules that the nodes fol- 
low to select where to send a certain packet. In all the 
topologies considered (ID, 2D and Cayley) the packets 
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follow paths of minimum length from their origin to their 
destination — open boundary conditions are set in both 
ID and 2D networks. In ID lattices and Cayley trees it 
is trivial to follow a minimum path because there is only 
one minimum path between two arbitrary nodes. In 2D 
lattices, however, there are many minimum paths. If a 
node can choose between two neighbors when sending a 
given packet, and both neighbors belong to a minimum 
path between the origin and the destination of the packet, 
one of them is chosen randomly with equal probability. 
This algorithm is indeed the simplest one and the inter- 
pretation of the results is clearer than for more complex 
routing algorithms. However, congested nodes will still 
have lower probability of receiving packets because of the 
definition of the quality of a channel. Therefore packets 
will avoid those nodes to some extent, as would happen 
in more sophisticated routing algorithms. 



III. RESULTS 

For certain parameters of the model (in particular 
£ = 1 jl6| ) and for similar models of traffic , a transi- 
tion from a free regime, where all the packets reach their 
destination, to a congested, regime where some packets 
are accumulated in the network, has been found. This 
transition is properly described by means of the follow- 
ing order parameter r](p) introduced previously pa: 



nip) 



t^oo pS At 



(4) 



In this equation AN = N(t + At) - N(t) and (...) in- 
dicates average over time windows of width At. These 
averages can be over one or many realizations, yielding 
the same result. Essentially, the order parameter repre- 
sents the ratio between undelivered and generated pack- 
ets calculated at long enough times such that AiV oc At. 
Thus, r\ is only a function of the probability of packet 
generation per node and time step, p. 

The power spectrum S(f) of the total number of pack- 
ets in the network N(t) = J2i n i 1S a ^ so used here to 
further understand the phase transition. By means of 
the power spectrum the behavior of the time correlations 
of the system can be studied. 

Let us study separately the critical case £ = 1 and the 
non critical cases £ < 1 and £ > 1. 



A. The critical case £ = 1 

A continuous transition between the free regime and 
the congested regime occurs in hierarchical networks for 
£ = 1 , as reported in Q . For small values of p, all the 
packets reach their destination and the total number of 
packets N(t) fluctuates around a finite value. In this case 
the order parameter is rj = 0. However, as p increases, a 
critical point p c is reached, where the fluctuations in N(t) 



become very large and the characteristic time of the sys- 
tem diverges (critical slowing down). Beyond this point, 
some packets remain undelivered and (N(t)) grows lin- 
early with t. The same qualitative behavior is observed 
here for ID and 2D lattices, although there are quanti- 
tative differences. 

From an engineering point of view it is interesting to 
study first the behavior of p c as the number of nodes in 
the network S changes, because it will provide valuable 
information about the scalability of the network. Note 
that N c = p c S is the maximum number of packets that 
the network can handle per time step and, thus, it is a 
measure of the capacity of the network. 

For the hierarchical Cayley tree, a mean field calcula- 
tion of p c was obtained in |p_6[] : 



z 3/2 1 



+ 1 (2-1)5' 



(•5) 



where the approximation holds for z m_1 3> 1. 

It is also possible to derive a mean field expression of 
p c for the ID lattice. Since the most congested node is, 
from symmetry arguments, the central one — the node at 
£ = L/2 — , the network will collapse when the amount 
of packets received by this central node is higher than 
the maximum number of packets that it is able to de- 
liver. Since in a large enough network it is safe to as- 
sume that the central node will be congested similarly to 
its neighbors, n^_i = nt = the maximum number 

of delivered packets should be 1. On the other hand, the 
number of packets arriving to the central node at each 
time step is the number of packets that are generated at 
each time step at the left half of the network and have 
their destination at the right half and conversely, this is 
pL/2. Then the critical condition is given by 



I = 



n 1D T 9 

P±± ^ pin = ± 
2 Pc L 



(6) 



These mean field expressions (|j) and (||) can be com- 
pared with simulations. Nevertheless, the fluctuations 
of N(t) become very large near p c and it is difficult to 
calculate the value of the order parameter. Instead, a 
susceptibility-like function x(p) can be defined by anal- 
ogy with equilibrium thermal critical phenomena, and 
used to estimate more accurately the value of the critical 
probability of packet generation, p c . The susceptibility x 
is related to the fluctuations of the order parameter by 



1 — >oo 



(7) 



where T is the width of a time window, and cr^(T) is 
the standard deviation of the order parameter estimated 
from the analysis of many different time windows of 
width T. Thus a calculation implies a long realization 
of N(t), its division into windows of width T, calcula- 
tion of the average value of the order parameter in each 
window and finally the determination of the standard de- 
viation of these values. The susceptibility shows (Fig. ||) 
a singularity at p c as T grows, as expected. 
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FIG. 2: Susceptibility for (a) a ID lattice with L = 100 and 
(b) a (5,4) Cayley tree, for different time windows T. The 
vertical dotted line corresponds to the mean field calculation 
of the critical point (equations (^) and respectively.) 

Fig. H shows that there is good agreement between ex- 
pressions (^|) and (|^) and the values of p c obtained nu- 
merically by means of susceptibility measures for differ- 
ent network sizes. 

For the 2D lattice it is more difficult to obtain even 
a mean field expression for p c . However, since for ID 
lattices and Cayley trees the scaling relation p c cx S^ 1 
holds — where S is the size of the network — , one may ex- 
pect the same behavior for the 2D lattice. Using the sus- 
ceptibility to numerically determine p c from simulations, 
one finds that this turns out to be incorrect. Although 
it is difficult to obtain a precise value of the exponent, 
Fig [| shows that it is close to 0.6 instead of 1.0: 



pl D oc 



s-°- 6 . 



(8) 



This result suggests that the existence of multiple paths 
to get from the origin to the destination has important 
consequences, not only in shifting the value of p c , but 
actually changing its critical scaling behavior. 

The behavior of the order parameter is studied next, 
ft is possible to derive an analytical expression for the 
simplest ID case where there are only two nodes that 
exchange packets. Since from symmetry considerations 
m = ri2, the average number of packets eliminated in one 
time step is 2, while the number of generated packets is 
2p. Thus p c = 1 and with the present formulation of 
the model it is not possible to reach the super-critical 
congested regime. However, p can be extended to be the 
average number of generated packets per node at each 
step (instead of a probability) and in this case it can 
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FIG. 3: Comparison between analytical (lines) and numeri- 
cal (symbols) values of p c obtained for the (a) ID and (b) 
the Cayley trees. The error bars of the numerical points are 
smaller than the size of the symbols. 



actually be as large as needed. As a result, the order 
parameter for the super-critical phase is i] — (p — l)/p. 
As observed in Fig. pL the general form 



V(p/Pc) = 



P/Pc ~ 1 

P/Pc 



(9) 



fits very accurately the behavior of the order parameter 
not only for this simple ID lattice with L = 2 or the Cay- 
ley tree , but also for any ID lattice. Two-dimensional 
lattices again deviate from this behavior, although the 
deviation is small. 

In particular, near p Cl equation (^) implies 



v(p) « (p - p c 



(10) 



and thus the critical exponent for the order parameter is 
equal to 1 at least for the ID lattice and the Cayley tree. 

The analysis of the power spectrum shows that in the 
sub-critical regime, i.e. in the free phase, the spectrum is 
well fitted by a Lorentzian characterized by a frequency, 
f Cl corresponding to exponentially decaying correlations 
with a characteristic time r, as depicted in Fig. ^. As 
p approaches p c , it is observed that f c — * and the 
power spectrum becomes l// 2 for the whole range of fre- 
quencies. Alternatively the characteristic time diverges 
as r — > oo (critical slowing down). This qualitative be- 
havior is common to all the topologies of the underlying 
network, as shown in Fig. o. 
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FIG. 4: (a) Susceptibility for a 2D 6 x 6 lattice, (b) Scaling 
of the critical probability of packet generation as a function 
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FIG. 5: Behavior of the order parameter in the case £ = 1. 
The solid line corresponds to the analytical calculation for two 
nodes exchanging information packets. Symbols correspond 
to simulations performed in ID, 2D and Cayley lattices. 



It is interesting to study how the characteristic fre- 
quency drops to for each network topology. Near the 
critical point, one expects the scaling behavior 



f c cx (p c - pY 



(11) 



The value of the critical exponent 7 can be estimated 
by fitting Eq. ([ll]) to values of f c (p) close enough to the 
critical point, as shown in Fig. ^. Note that we fit p c and 
7 simultaneously. This procedure yields very accurate 
values of p c but the values of 7 are subject to large errors. 
Fig. H yields 7 « 0.9 for a ID network with L = 100, 
7 w 2.5 for a 2D network with L = 7, and 7 w 2 for a 
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FIG. 6: Left: Log-Log plot of the power spectrum of N(t) 
for different topologies and different values of the control pa- 
rameter e = (p c — p)/p c for the ID case (L = 100), 2D case 
(L — 6), and Cayley (7, 5). Power spectra have been obtained 
averaging over 100 realizations of N(t). Dotted lines repre- 
sent a power law with exponent -2. Right: Characteristic 
frequency as a function of the control parameter e for the dif- 
ferent topologies. As observed, the characteristic frequency 
tends to as p — > p c following a power law. The straight lines 
correspond to fittings of equation (|ll|). 



(7,5) Cayley tree. 

The determination of 7 is interesting not only from 
an academic point of view, but also from an engineering 
perspective. Indeed, this exponent is related to diver- 
gences of other relevant quantities near the critical point. 
Any characteristic time r — the average time to deliver a 
packet, for instance — will diverge as 



T OC (p c - p) 7 

and similarly the total number of packets 

cx (p c -py 



since 



— = pS 

T 



(12) 



(13) 



(14) 



from Little's law This law states that, in steady 

state, the number of delivered packets and the number 
of generated packets are equal. 
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The estimation of 7 is particularly interesting in elec- 
tronic communication protocols. Indeed, Eq. ( |l2[ ) is used 
to determine the waiting time before a packet is consid- 
ered lost in the network and therefore sent again |l3| . In 
practice, the exponent 7 = 1 predicted by classical queue 
theory is assumed, while our current results suggest 
that more complex settings can lead to exponents even 
larger than 2. 

B. Non critical cases £ < 1 and £ > 1 

We have shown that the number of packets delivered by 
node i is n\~ and thus, when £ < 1, it increases with the 
number of packets that this node accumulates. It is dif- 
ficult to imagine a real scenario with this characteristic. 
However, this case has been included to understand the 
critical behavior when £ = 1, i.e., to show the relation- 
ship between criticality and the amount of packets that 
can be delivered when load increases. As a consequence 
of the increase of the deliver capability with the load, the 
transition to collapse will never occur because, at some 
point in time, the number of accumulated packets will 
be large enough and the number of delivered and created 
packets will balance each other. Thus, the order param- 
eter will be zero for any value of the control parameter p, 
and the correlations will decay exponentially. As shown 
in Fig. |?], the characteristic frequency tends asymptoti- 
cally to /* as p increases. This asymptotic value depends 
on the size of the system. 

For a ID lattice with a high density of packets (p — > 1), 
the number of packets that are delivered by a node is 
n\~^ while the number of packets that are being delivered 
to this node is proportional to L (for instance, for the 
central node, this number is simply pL/2). Therefore, 
Hi oc L^/C 1 -?). The total number of packets is N = 
Yli n i ~ L 1+1 /(l-£) and according to Little's Law 

/ c * ex ^ oc (15) 

On the other hand, for p — > the scaling of the char- 
acteristic frequency is given by 

f c « L" 1 (16) 

since the packets success to jump from one node to the 
next at each time step, and therefore the characteristic 
time is directly the average path length between nodes. 

Therefore, although there is no phase transition in this 
case £ < 1, there is a cross-over from a low density be- 
havior to a high density behavior, as shown in Fig. |t]. 
This crossover is also observed in 2D lattices and Cay ley 
trees. 

The phase transition observed for £ = 1 is recovered 
when £ > 1. Above a certain threshold, some packets 
are accumulated in the network and the order parameter 
differs from 0. However, the number of packets delivered 




Control parameter, p 

FIG. 7: Characteristic frequency f c as a function of the prob- 
ability of packet generation p, for £ = 0.2 and different sizes 
of a ID lattice. As observed, f c never becomes as happens 
in the critical £ = 1 case. Inset: Characteristic frequency 
at p — *• 0, f° (squares), and characteristic frequency at large 
p, /* (circles). The lines represent the fittings provided by 
equation (jl6|) /° oc L~ l , and equation ( |l5| ) /* oc L^ 1 ^ 1 " 5 ', 
respectively. 




FIG. 8: Congestion nuclei formation for large 2D lattices of 
L = 200, in the non critical case £ > 1. Dark regions represent 
regions with small congestion levels while bright regions corre- 
spond to highly congested regions. Left: £ = 5 and p = 0.001. 
Right: £ = 2 and p = 0.01. 



by a node i in this case £ > 1 decreases with the number 
of packets accumulated at that node. Therefore, when 
some packets are accumulated, n, grows and finally no 
packets are delivered at all. Thus suddenly above the 
transition, which is discontinuous, the order parameter 
becomes 1. 

The change in the order of the phase transition affects 
the spreading of the collapse over the network. In the 
critical case f = 1, the collapse starts at the most central 
node and then spreads from this point to the rest of the 
network. In this case £ > 1, the reinforcement effect — 
the fact that the more collapsed a node is, the more col- 
lapsed will get in the future — leads to the formation of 
many congestion nuclei generated by fluctuations, that 
spread over the whole network. Figure ^| illustrates the 
formation of these congestion nuclei for 2D lattices with 
£ = 5 and p = 0.001, and £ = 2 and p = 0.01 
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IV. DISCUSSION AND CONCLUSIONS 

A collection of models recently proposed for hierarchi- 
cal networks ]l6| has been examined in detail for several 
network topologies including ID and 2D lattices, to char- 
acterize the phase transition to collapse. New congestion 
scenarios have been analyzed. 

First, it has been shown that the congestion behavior 
is governed by the ability of agents to deliver information 
packets when their load increases. When agents deliver 
packets at a constant ratio independently of their load — 
number of packets to deliver — , a continuous transition as 
the one reported in jl6) for Cayley trees is also observed 
for ID and 2D lattices, the order parameter being the 
fraction of accumulated packets. When the number of 
delivered packets increases with the load, no phase tran- 
sition is observed. Conversely, when the number of deliv- 
ered packets decreases with the load, the order parameter 
jumps from zero to one and the transition becomes dis- 
continuous. These different behaviors are tuned by a sin- 
gle parameter £, which is the exponent that determines 
how the capability of nodes evolves with the number of 
accumulated packets. When £ < 1 there is no transi- 
tion to the congested regime, for £ = 1 the transition is 
continuous and for £ > 1 the transition is discontinuous. 
Note that the continuous transition (reported in (l6) and 
in some models of queues [P7| ) is only a particular 
case between a no-congestion behavior and a discontin- 
uous transition behavior. Thus, the critical behavior is 
intimately related to the independence between load and 
deliver capability of nodes. 

These different behaviors have been analyzed sepa- 
rately. The critical case £ = 1 presents the most interest- 
ing behavior. To properly understand the physics of the 
collapse process it is necessary to know how the system 
approaches the critical point. The order parameter be- 
comes useless in this region due to the emergence of large 
fluctuations, and then we have defined a susceptibility- 
like function that is related, by analogy to equilibrium 
critical phenomena, to fluctuations of the order parame- 
ter. This quantity shows a peak at the critical point and 
allows to study the dependence on the size and topology 
of the system. For ID lattices and Cayley trees, it has 
been shown (by means of simulation and also with mean 
field calculations) that p c scales as the inverse of the size 
of the system, p c cx S . However, for 2D lattices a dif- 



ferent scaling p c cx S~ 6 is obtained. It is suggested that 
the existence of multiple paths from the origin to the des- 
tination of the packets is the responsible for this change 
in the scaling behavior. On the other hand, the study of 
the power spectrum of the number of the packets in the 
network as a function of time has provided valuable in- 
formation about the variation of temporal correlations in 
the system. In particular, below p c correlations decay ex- 
ponentially, with a characteristic time r that diverges at 
p c . The exponent of this divergence is significantly larger 
than what is expected from queue theory. This can be 
relevant to accurately forecast the temporal behavior in 
complex communication settings. 

For £ < 1, no phase transition is observed. Instead, 
a cross-over from a low-density to a high-density regime 
occurs. In the low density regime, the characteristic fre- 
quency is determined by the average distance between 
nodes, while in the high density regime the characteristic 
frequency is determined by the capability of agents to de- 
liver information packets. In ID lattices, the cross-over 
changes the scaling from f c cx L^ 1 at low densities to 
f c oc L^ 1 ^ 1 ^^) a t high densities. 

Finally, when £ > 1 the transition to the congested 
regime is discontinuous, since the number of delivered 
packets decreases with the load and for long times no 
packets are delivered at all. Thus the order parameter 
jumps from to 1. Moreover, due to the reinforcement 
effect — the fact that the more collapsed a node is, the 
more collapsed it will get in the future — leads to the for- 
mation of many congestion nuclei generated by fluctua- 
tions, that spread over the whole network. The existence 
and distribution of such congestion nuclei can also be of 
interest from an engineering point of view. 
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